Electroweak interactions and the muon g — 2: bosonic two-loop effects 
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We present a detailed evaluation of the bosonic two-loop electroweak corrections to the muon's 
anomalous magnetic moment. We study the Higgs mass dependence and find agreement with a 
previous evaluation in the large Higgs mass limit. We find a^ w bos (two- loop) = (— 22.2±1.6) x 10 -11 , 
for 114 GeV < Af m gg S < 700 GeV. 
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I. INTRODUCTION 



The anomalous magnetic moment of the muon, a M = (<7„ — 2) /2 has recently been determined 
with a very high precision. A series of measurements with positive and negative muons by the E821 
Collaboration at the Brookhaven National Laboratory resulted in the present average value 



(116 592 080 ±60) x 10" 



(1) 



At this level of a one half part per million precision, this quantity is sensitive to subtle effects predicted 
by the Standard Model (SM), including five-loop quantum electrodynamics (QED), hadronic vacuum 
polarization and light-by-light scattering, and electroweak interactions at two-loops. 

Obviously, — 2 may be affected also by interactions beyond the SM. For this reason, many 
researchers have analysed g^ — 2 in various models of new physics and performed sophisticated studies 
of the SM effects which are an irreducible background in the search of unknown phenomena. A 
summary of the SM prediction can be found, for example, in the recent studies and reviews 0, 0, 0, • 

At present, the experimental result in Eq. (JIJ exceeds the SM prediction by about 2.6 times the 
combined theoretical and experimental uncertainty. This tantalizing discrepancy may be due to an 
effect of a new interaction, perhaps supersymmetry. However, it is important to scrutinize the SM 
prediction before a conclusion can be made. The present paper is devoted to a reevaluation of the 
largest part of two-loop electroweak diagrams, namely those without closed fermionic loops. 

Electroweak one-loop corrections to g^ — 2 were among the first quantum effects studied in the 
renormalizable electroweak theory, in 1972 They were found to increase g M — 2 by 



EW 
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1 - 4 sin 2 



\ M 2 



195 x 10~ n , 



(2) 



M w jM\. The large mass 



where G M = 1.16639(1) x 10~ 5 GeV^~ 2 and weak mixing angle sin 2 9 W = 1 
parameter M represents the mass of a W, Z, or a Higgs boson. 

This effect was too small to be measured in the then ongoing CERN experiment. The desire to 
observe it motivated in part the latest Brookhaven effort. 

Two decades after the first electroweak result, Kukhto et al. 0] found that an additional virtual 
photon may significantly modify the one-loop value in Eq. ©. They estimated that effect as a -22% 
reduction - surprisingly large for hard virtual photons. This and an analogous reduction of the rare 
muon decay fi — > &y [l2j are due to the large anomalous dimension of the dipole operators such as 
fia^nF^ and eo^pF^. 

That large effect found in a subset of two-loop electroweak contributions was similar in size to 
the design precision of the E821 experiment. It therefore appeared as warranted, even necessary, 
to evaluate the complete two-loop result - a calculation that had not been performed before in the 
electroweak theory for any other observable. 

In 1995 a complete set of 1678 electoweak two-loop diagrams for g^ — 2 was generated by a computer 
system [13J . However, such large number of diagrams, many of which are divergent, could not be 
numerically calculated, at least at that time. Fortunately, it turned out that the majority of those 
diagrams are strongly suppressed by extra factors of the muon-to-intermediate-boson mass ratio and 
can be neglected. Thus, a complete two- loop result was found 0|- The numerical value was found to 
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be dominated by large logarithms arising due to photon exchanges and was, somewhat accidentally, 
close to the value found in its first studies 0, 0] . 

In later work, the renormalization group equation was employed to estimate higher-order logarithmic 
effects, which however are not sizable [l6l Il7| . 

The result of |14| was obtained in an approximation assuming that the Higgs boson is sufficiently 
heavier than the W and Z bosons such that M\y,z /Mjj can be used as an expansion parameter. A 
recent study |l8j | relaxed that approximation and found the value of two-loop contributions valid also 
for a light Higgs (moreover, the two- loop supersymmetric effects were evaluated in that paper). 

In this paper, we reanalyse the two-loop electroweak effects. Our goal is to check the previous 
results and study the Higgs mass dependence. We present our results in a semi-analytical form. That 
is, the dependence on the Higgs mass is presented analytically, while some parts of expressions that 
depend only on the well-known particle masses are, for the sake of brevity, evaluated numerically. A 
method employed for obtaining analytical results is also described in some detail. 

In this work we focus on diagrams with only bosonic loops (no closed fermion loops). The fermionic 
subset of corrections was studied separately |l9ll20l l2ll | . It was recently subject of an interesting the- 
oretical controversy which seems to be settled now (see P3 for a thorough discussion and references). 

In the next Section we briefly explain the asymptotic operation which is the main technical tool 
used to obtain the analytical result for g^ — 2. Section IlIII presents partial results for various groups 
of contributing diagrams. First, we divide all diagrams into five subsets of topologically equivalent 
diagrams. Then we give analytical or semi-analytical result for each topology. Finally we discuss 
the renormalization procedure and evaluate necessary counterterms, provide the final numerical result 
and compare it with the results of Ref. 0, . 



II. TECHNIQUES: ASYMPTOTIC OPERATION AND RESULTING INTEGRALS 

The method of asymptotic operation (see for a review and references) has been an invaluable 
tool in numerous recent studies of effects involving various energy scales. In the present problem, the 
mass of the muon m sets the soft scale, and the masses of the W, Z, and Higgs bosons - the hard 
scale: Mw,z.h 3> m. As we explained in the Introduction, the original study 14) assumed in addition 
Mh 3> Mw,z- Here, we will not assume any hierarchy between Mw,z and Mjj at the price of certain 
complication of our results. 

We would like to explain here the basic principles of asymptotic operation. Instead of attempting 
a rigorous derivation or even a rigorous exposition, we will use a simple example to elucidate the 
method. 

In this study, we are interested in two-loop Feynman graphs G, which have a soft-scale (~ m) 
external momentum and involve internal lines with both soft-scale and hard-scale (~ M) virtual 
particles. The exact value of the two loop integrations is a (possibly very complicated) function of 
m/M . However, for our purposes it is entirely sufficient to know only an expansion of that function up 
to m 2 /M 2 . The purpose of the asymptotic operation is to obtain the desired order of that expansion 
without having to compute the whole function. 

For our purposes, the action of the asymptotic operation on a Feynman graph G may be described 
with the following formula, 

As o G = V {roh)x [G\h] . (3) 

^ — ' h 

Let us first of all explain the notation. As o G is an expansion of the exact value of G in powers and 
logarithms of m/M . h are the subgraphs of G in which all loop momenta are considered as hard. 

Instead of giving an exact definition of which h are relevant, we use as an example the diagram shown 
in Fig. ^ In that Figure, G is shown in picture (a). Subgraphs h are shown in (b)-(e), and range 
from a single heavy particle line to the whole graph, r is the Taylor expansion operator, expanding a 
subgraph in all soft-scale parameters such as the light-particle masses and external momenta of that 
subgraph (the latter include loop momenta of those loops that are not part of the subgraph) . Finally, 
G\h denotes the graph G with the subgraph h contracted into a point. The physical intuition for that 
concept is such that all masses and loop momenta in G\h are soft; thus, from their point of view, the 
hard-scale subprocess in h occurs at a very short-distance scale and can be contracted to a point. In 
other words, the subgraphs h create effective vertices. 

What do we gain from this construction? We completely separate soft and hard scales for the 
purposes of the integration. Thus, the resulting integrals cannot contain both m and M - they 
factorize into a part that depends only on the soft scales and another that depends on the hard ones. 



y 

(a) (d) (e) 

Figure 1: An example of applying the asymptotic operation to a graph G, shown in (a), (b)-(e) show the full 
set of subgraphs, as described in the text. Propagators through which hard momenta flow are indicated with 
solid lines - they constitute a subgraph h. Dashed lines contain only soft-scale momenta and masses - they 
constitute G\h. 



No non-trivial functional dependence on m/M can arise and we obtain the expansion to a given order 
in m/M by taking a sufficient number of terms in the Taylor expansion of the subgraphs. 

The asymptotic operation greatly simplifies the types of integrals that we need to evaluate, as shown 
in Fig. @. One- loop on-shell integrals and one-loop massive integrals, pictured in Fig. are 
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Figure 2: Basic types of loop integrals needed for the evaluation of two-loop electroweak corrections to <? M — 2: 
(a) One-loop light-mass on-shell integral multiplied by a one-loop vacuum integral with a large mass; (b) 
Two- loop vacuum integral (some propagators may be massless (M — 0)); (c) Two-loop light-mass on-shell 
integral. 



trivial. 

Two-loop integrals in Fig. J5J(b), dependent on arbitrary masses, were computed in Ref. |23t 
l24j . By means of integration by parts identities |25| the initial integral with arbitrary powers of 
propagators can be reduced to the integral with all propagators in power one. For the case when all 
three propagators are massive, we encounter only diagrams with two different masses present. In this 
simplified case, the general result for such integral, in the dimensional regularization, becomes 
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where d = 4 — 2e, z 
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and the function $ is defined as 
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d01n|2sin(0/2)| 



-4Li 2 (0 + 21n 2 £ - ln 2 4z + — 



(4) 



(5) 



(6) 



for z > 1. 

For the case when one propagator is massless we introduce the function I Rn , which is the finite part 
of the integral I 



M c M a M b 



for the case of M c — 0; denoting R — -r£ we have 



1 



I nn {R) = -^(l + R)+3RlaR- (1 -H)Li 2 (l - R) - — In 2 R. 



(7) 
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Finally, two-loop on-shell integrals shown in Fig. J2J( C ) can be found, for example, in Ref. p6| . 

To calculate the renormalization counter terms we shall also need one-loop massive on-shell self- 
energy integrals shown in Fig. J3J). They are expressed as linear combinations of integrals of the type 



Mi 



M2 



Figure 3: One-loop massive on-shell self energy integral 



(we use p 2 = — M 2 ): 

S ab (M!,M 2 , M) = ^ ~^ ± £ \ f dx 77— 5— y . (8) 

r(6)r(l + e)7 [(1 - X )(M 2 - xM 2 ) + xM 2 f- 2+£) 

These integrals can be represented as series in e and can be easily computed to the necessary order. 
We introduce special functions S^(M\, M 2 , M) which represent the finite parts of ©. We shall need 
only the following cases, 

S£l(M 1 ,M 2) M) = f dx{ [(I - x) (M 2 - xM 2 ) + xM 2 ] [in ((1 - x) (M 2 - xM 2 ) + xM 2 ) - l] } , 
Jo 

Sil{M 1 ,M 2 ,M) = - I dxx a In [(1 - x) (M 2 - xM 2 ) + xM$] , (9) 
Jo 

with a taking one the values 0, 1 or 2. 



III. BOSONIC ELECTROWEAK TWO-LOOP CONTIBUTIONS IN THE SM 



In this section we apply the technique of asymptotic operation, as described in Sec. [nj to evaluate 
all two-loop bosonic contributions. We divide the diagrams into subsets of various topologies and 
mass assigments, and present detailed results for each subset. We display the dependence on the 
Higgs mass in the analytical form. We also provide details of the renormalization procedure. 

Following [3| j we use the 't Hooft-Feynman non- linear gauge • We choose it in such way that 
we eliminate the vertex jW^G^. Since — 2 involves an external photon in every diagram, such 
choice greatly reduces the number of diagrams. The set of diagrams we have to consider contains all 
two-loop diagrams one can compose in SM with the exclusion of pure QED diagrams and diagrams 
with closed fermion loop; we drop all diagrams with more than one scalar coupling to the muon line, 
since each such coupling introduces an extra factor jj^- Taking advantage of the mirror symmetry 
we reduce the number of diagrams that must be individually calculated to 138; in addition, some 
one-loop diagrams must be evaluated for the renormalization. 



A. Two-loop topologies and their evaluation 



We divide all two- loop diagrams into five topological types, as shown in Fig. Possible insertions 
of the external photon vertex are indicated with a circle cross, and wavy lines stand for either scalar 
or vector bosons. Each subset groups diagrams with similar properties with respect to the asymptotic 
operation and thus can be calculated using the same algebraic code (we use FORM |28| for most 
algebraic operations). 

Applying the general formula J3J for the asymptotic operation and computing all resulting integrals, 
we obtain the following results for the finite parts of various groups of diagrams: 

T1a,b,c + T2a,d,e + T3a,b,c,d,g + T4a,d,e,j + T5 
a 2 

= |680c 4 - 362c 2 - 363 + 6 ([27 - 30c 2 l mm 2 + [54c 2 - 56c 4 l lnM 2 ^ 
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Figure 4: Diagram topologies and assignments of masses to their lines. 



[56c 4 - 84c 2 + 27] InMf)} - a 2 ^- (4.6 - O.lQTIn-^-) 
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In the above expressions we have dropped the divergences (which cancel in the final sum with the 
counterterms) . We also express all masses in units of // = 1 GeV. The floating point coefficients in 
HI U| l were rounded to provide the precision of 10 -11 of the final result. They were obtained using 
M w = 80.423 GeV and M z = 91.1876 GeV. The functions $(x) and I Rn {x) are defined by means of 
(£0 El and 01 respectively. We also use the notation A=|^A=ifa,c=ir^ands = Vl - c 2 . 

The results for topologies T1 A , T1 B , Tl c , T2 A , T2 D , T2 E , T3 A , T3 B , TZ C iTZ D , T3 G , T4 A , TA D , 
T4e, T4j, T5, were collected into a single formula because the corresponding diagrams do not contain 
Higgs propagators and thus can be evaluated without any assumptions about Mjj- 



B. Renormalization counterterms 



Conceptually, the renormalization procedure for our calculation is identical to the one we used in 
Ref. [l4(- The only technical difference is that in our calculation we computed one-loop massive 
self-energy diagrams precisely, rather than in an expansion to the order of sin 6 6w- In this section we 
briefly review our renormalization procedure and list specific expressions for the W and Z self-energies 
and for other renormalization constants. 

The counter terms for the two-loop EW corrections are generated by renormalizing vertices and 
propagators in one-loop diagrams. In the non-linear 't Hooft-Feynman gauge there are just three such 
diagrams: the Schwinger's QED diagram with a photon loop, its analog with the photon replaced by 
a Z boson, and a diagram with two W bosons. 
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The result for Schwinger's diagram in dimensional regularization is given by 



Schw 



= ii[l +£ (4-W)] 



(19) 



For the calculation of the counterterms we also need the contribution of that diagram with one internal 
muon propagator squared. It is 
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Using (|19fl and (|2U|) . and bearing in mind that since for the calculation of the EW corrections we 
droped the two-loop photonic corrections to the Schwinger diagram, and thus have to subtract the 
photonic contributions from muon mass and wave function renormalization constants as well, we find 
the counterterm to the a/ at the two-loop level, 
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The corresponding diagram with one muon propagator squared gives 



M% 4tt 



(g 2 v - 9\) 



1 



— e 



lnm 2 



(23) 



Combining expressions 1|22|) and (|23fl we obtain the counterterm generated by the Z-loop diagram, 
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Finally, the one-loop contribution of the diagram with two W bosons is 
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Using this expression we obtain a relatively simple counter term because there is no internal muon 
line and there are only left-handed couplings, 
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The constants appearing in the expressions (|19(I - (|2G|I are defined as follows: 
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and 5m, SZff' are the muon mass and wave function renormalization defined as 
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and SM W and <5M§ are mass renormalization constants. 

The expressions (|21l 1241 126|l provide all counter terms necessary to renormalize two-loop bosonic 
corrections to a 
terms, 

„,2 



p. Taking the sum of the finite parts of these expressions we obtain for the counter 
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where, as in expression l|10[l . the floating point coefficients were rounded to provide the precision of 
10~ n of the result and do not depend on Mw and M z chosen within experimentally accepted interval. 
All special functions S a l ™ are defined in Eq. (JSJ. 
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C. Two- loop corrections to g — 2 

Our result (not yet in its final form - see the next section) for the two-loop bosonic contribution to 
gy, — 2 is obtained by adding up all diagrams, Eqs. |TDJ)— JTSJl, and counterterms i|30|l . 

a™ bos (two - loop) = Tl A , B ,c + T2 AtD<E + T3 a ,b,c,d,g + TU,d,e,j + T5 + T2 B 

+T2 C + T2 F + T3e,f + TA B + T4 C + Tip + T4 G + CT. (31) 

Its numerical value, normalized to the one-loop correction giving in Eq. ©, is plotted in Fig. JSJ) 
with a solid line, for a range of the Higgs boson mass 50 GeV < Mjj < 700 GeV. Also plotted, with a 
dashed line, is the approximate result of Ref. [l4j. We see that both results coincide well for the Higgs 
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Figure 5: 
Ref. fijl. 



(two — loop) 



as a function of Mh, expressed in percents. Dashed line represents the result from 



EW / i 
a l± (one — loop) 

solid line is the result of this work. The vertical dotted line shows the lower limit for the Higgs 
boson mass from direct searches, 114 GeV. 



boson heavier than about 200 GeV. The only difference in the large Higgs mass region is due to an 
additional approximation made in Ref. , where the difference of the W and Z masses was treated 
as a small perturbation (first four terms in an expansion in sin 2 6\v were retained there) . Above 
Mh = 250 GeV, the relative difference between two results is less than 0.1%, well within the precision 
of the result in Ref. 14j. However, the results differ strongly in the low Higgs mass region. The result 
of Ref. 14j is not valid in that re gion since it was obtained under the assumption M\y,z "C Mh- We 
see that while the result of Ref. [Tj| seems to grow strongly for the light Mh, the actual bosonic 
correction (solid line) remains moderate. 



D. Reparametrization in terms of G M 

Eq. (|31fl gives a finite result for the two-loop correction, expressed in terms of the fine structure 
constant a, the W mass, and the weak mixing angle. Some of the corrections computed in this way 
are universal for all weak processes and it is convenient to include them in the lower order result by 
expressing it in terms of the Fermi constant . This amounts to the substitution, to be made in all 
components of Eq. iloT|) , 
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with 



A, _„«,_,«!_ + « (. + I=£| h< ?. ,33, 

s M^, 4tt s 2 V 2s 2 ' 



This transformation decreases the central value of the ratio of the two-loop to one-loop corrections by 
up to 7% (depending on the Higgs mass) and, in addition, reduces its uncertainty, since G M has been 
measured to much better precision than Mw and s. 



IV. CONCLUSIONS 



We have presented a detailed evaluation of the two-loop bosonic corrections to — 2. Our result 
confirms the previous approximate evaluation in the limit of the heavy Higgs boson |14j . Our final 
number for this correction is 

a™ bos (two - loop) = (-22.2 ± 1.6) x 10~ n , (34) 

where the central value was computed for the Higgs mass of 200 GeV, and the error encompasses 
the interval 50 GeV < Mh < 700 GeV. This is in agreement with the numerical results of an exact 
study in Rcf. 18] . In addition to the numerical result we have provided here a number of semi-analytic 
intermediate results. They give insight into details of our calculation, allow future checks, and simplify 
the evaluation of the correction for any Higgs boson mass. 

The uncertainty in Eq. I|34(l is due to the unknown Higgs mass and is slightly reduced in comparison 
with previous estimates based on the approximate evalulation. The main variation of the correction 
occurs in the relatively low Higgs mass range. Thus, when more stringent limits on My are obtained, 
that uncertainty will quickly decrease. The main uncertainty of the electroweak contributions will 
then be due to electroweak hadronic effects, as discussed in [l6j . 



Acknowledgments 

This research was supported by the Natural Sciences and Engineering Research Council of Canada. 
TG gratefully acknowledges support by the Golden Bell Jar Department of Physics Graduate Schol- 
arship. 



[1] G. W. Bennett el ai, Phys. Rev. Lett. 92, 161802 (2004). 

[2] M. Passera, J. Phys. G31, R75 (2005). 

[3] K. Melnikov and A. Vainshtein, (2003), hep-ph/0312226 

[4] M. Davier and W. J. Marciano, Ann. Rev. Nucl. Part. Sci. 54, 115 (2004). 

[5] A. Czarnecki, Nucl. Phys. B (Proc. Suppl.) 144, 201 (2005). 

[6] K. Fujikawa, B. W. Lee, and A. I. Sanda, Phys. Rev. D6, 2923 (1972). 

[7] R. Jackiw and S. Weinberg, Phys. Rev. D5, 2396 (1972). 

[8] G. Altarelli, N. Cabibbo, and L. Maiani, Phys. Lett. B40, 415 (1972). 

[9] I. Bars and M. Yoshimura, Phys. Rev. D6, 374 (1972). 
[10] W. A. Bardeen, R. Gastmans, and B. E. Lautrup, Nucl. Phys. B46, 315 (1972). 
[11] T. V. Kukhto, E. A. Kuraev, A. Schiller, and Z. K. Silagadze, Nucl. Phys. B371, 567 (1992). 
[12] A. Czarnecki and E. Jankowski, Phys. Rev. D65, 113004 (2002). 
[13] T. Kaneko and N. Nakazawa, (1995), hep-ph/9505278 

[14] A. Czarnecki, B. Krause, and W. Marciano, Phys. Rev. Lett. 76, 3267 (1996). 

[15] E. A. Kuraev, T. V. Kukhto, and A. Schiller, Sov. J. Nucl. Phys. 51, 1031 (1990). 

[16] A. Czarnecki, W. J. Marciano, and A. Vainshtein, Phys. Rev. D67, 073006 (2003). 

[17] G. Degrassi and G. F. Giudice, Phys. Rev. D58, 053007 (1998). 

[18] S. Heinemeyer, D. Stockinger, and G. Weiglein, Nucl. Phys. B699, 103 (2004). 

[19] A. Czarnecki, B. Krause, and W. Marciano, Phys. Rev. D52, 2619 (1995). 

[20] S. Peris, M. Perrottet, and E. de Rafael, Phys. Lett. B355, 523 (1995). 

[21] M. Knecht, S. Peris, M. Perrottet, and E. de Rafael, JHEP 11, 003 (2002). 

[22] F. V. Tkachev, Sov. J. Part. Nucl. 25, 649 (1994) hep-ph/9701272 . 

[23] A. I. Davydychev and J. B. Tausk, Nucl. Phys. B397, 123 (1993). 

[24] F. A. Berends and J. B. Tausk, Nucl. Phys. 421, 456 (1994). 



11 



[25] F. V. Tkachov, Phys. Lett. BlOO, 65 (1981). 

[26] J. Fleischer and M. Yu. Kalmykov, Comput. Phys. Commun. 128, 531 (2000). 

[27] K. Fujikawa, Phys. Re v. D7, 393 (1973) 

[28] J. A. M. Vermaseren, math-ph/0010025 



